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POSITIVE  DEFINITE  SPECTRAL  ESTIMATE  AND  STABLE 
CORRELATION  RECURSION  FOR  MULTIVARIATE  LINEAR 
PREDICTIVE  SPECTRAL  ANALYSIS 


INTRODUCTION 


A generalization  of  burg's  algorithm  for  spectral  analysis  to  the 
multivariate  case  was  the  subject  of  an  earlier  report.^  All  the  desir- 
able properties  of  the  univariate  case  were  shown  to  hold  true,  except 
that  it  was  not  proven  that  the  residual  matrix  was  positive  definite, 
nor  that  the  correlation  recursion  was  stable.  Both  of  these  assump- 
tions can  be  affirmed  by  drawing  on  the  results  in  Strand^  and  Burg.^ 

In  addition  to  affirming  these  two  assumptions,  this  report  contains 
a modified  and  updated  FORTRAN  program  that  supersedes  the  program  pre- 
viously reported.!  The  modified  program  incorporates  some  more-explana- 
tory format  statements,  the  calculation  of  the  (normalized)  correlation 
matrices  via  recursion,  and  the  aliased  (normalized)  correlation  matrices 
by  means  of  a Fast  Fourier  Transform  (FFT) . 

This  report  is  a sequel  to  an  earlier  report.!  In  order  to  elimi- 
nate duplication,  that  report  is  referenced  for  background  information, 
a list  of  symbols  used,  and  processing  technique.  We  shall  draw  freely 
on  that  report;  for  example,  equation  (5)  of  the  earlier  report  will  be 
denoted  by  (5) . ! 

POSITIVE  DEFINITE  RESIDUAL  MATRIX 

The  (p-l)-th  order  forward  residual  matrix,  Up_j,  was  defined  in 
equation  (95).!  ^ish  to  show  that  Up  is  positive  definite;  the 
following  proof  is  based  on  reference  2,  equations  (3.25-3.32). 

From  equation  (H-5),!  we  have,  using  the  Hermitian  property  of 
Up  and  Vp, 

^ V^,  ; (1) 

and  from  equation  (137),!  eliminating  B^P^^, 


(2) 
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Notice  that  we  have  made  specific  use  of  the  inverse  weighting  in 
equation  (136). ^ Substituting  equation  (2)  into  equation  (1),  we  find 


therefore, 


Taking  the  conjugate  transpose  of  both  sides  of  equation  (4)  and 
using  equations  (106)1  and  (114)1  yields 

= s;!f -2^;:  . 

Adding  equations  (4)  and  (5)  together  and  multiplying  by  -1,  there 
follows  ^ 

= -2  \5^-  A’:’  5^"-  s;*?  A»’\  At  s;f  At"] =-  2Ep  > 

the  last  identity  was  derived  from  equation  (113).^ 


Define 


Mf.,  - 4::  5;?r . 


Then  equation  (6)  becomes  simply 

Mfl.  4 4 M,.,  = - 2 Ef . (8) 

Now,  Ep  is  Hermitian  and  positive  definite  (see  equation  (112)^); 

also,  is  Hermitian  and  positive  definite  (see  equation  (114A)^). 

We  assume  that  U ^ is  positive  definite.  Then,  Uplj  is  positive 

definite,  and  so  U‘l  must  have  all  its  eigenvalues  positive 

p-1  p-1 

WToTTKe-positive  definite  statements  should  be  qualified  with 

the  proviso  "with  probability  1." 
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(see  appendix  A).  As  a result,  Mp_l  has  all  its  eigenvalues  negative, 
making  it  a stable  matrix  (reference  4,  page  270).  Therefore,  the 
solution  of  equation  (8)  exists  and  is  unique  (reference  5,  equation  3). 

According  to  reference  4,  page  278,  problem  3,  there  exists  a posi- 
tive definite  solution  of  equation  (8)  for  Up.  Therefore,  there  is  a 
unique  positive  definite  solution  of  equation  (8)  for  U . Since 


H 


(9) 


(from  equations  (95)^  and  (82)^)  is  positive  definite,  the  assumption 
above,  that  Up-i  is  positive  definite,  can  be  justified  by  induction. 

In  summary,  the  residual  matrix  Up,  calculated  by  means  of  equa- 
tion (105)i  or  (181), i is  positive  definite.  The  quantity  Vp  is  also 
positive  definite;  the  equation  analogous  to  equation  (6)  is 


(10) 


and  all  the  comments  above  apply  directly.  It  is  worth  repeating  that 
the  positive  definite  conclusion  on  Up  and  Vp  holds  for  the  specific 
inverse  weighting  indicated  in  equation  (136)1;  whether  it  also  holds 
for  other  weightings  is  unknown. 


STABLE  CORRELATION  RECURSION 

The  correlation  recursion  is  given  in  equation  (164)^  according  to 


where  superscript  p has  been  added  to  the  correlation  matrices  to  indi- 
cate specifically  their  dependence  on  the  p-th  order  predictive  filter; 
and  starting  values  have  been  defined,  as  in  equation  (D-3),l  namely. 
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= 'R-  . HI  * r- 


(12) 


The  latter  quantities  in  equation  (12)  are,  according  to  equation  (78A),1 
solutions  of 


ll»  I 


, Isoisp 


(13) 


Combining  equations  (11)  through  (13),  we  have 


(14) 


We  will  show  that  recursion  (11)  is  stable;  that  is,  we  will  show 
that  (the  elements  of)  matrix  does  not  tend  to  infinity  as  m tends 

to  infinity,  with  p fixed.  The  proof  is  an  extension  of  reference  3, 
section  I11.C.2  (which  was  for  known  correlation),  to  fit  the  unknown 
correlation  case. 

We  have,  from  equations  (82)1  (gOA),^  respectively, 

v” 


-for 


(15) 


For  a given  value  of  p,  define  the  (m  + 1)  x (m  + 1)  block  Toeplitz 
matrix 


(R 


ip) 


ic  T»r 


■R'f' 

•Hi 


(16) 
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If  m ^ p,  the  entries  in  equation  (16)  are  according  to  equation 
(12),  whereas  if  m > p,  the  entries  are  those  generated  by  equation  (11). 
It  follows  immediately,  from  equations  (16)  and  (12),  that 

OHa  = aC  if  s r- 

Hie  s,t-th  block  of  in  equation  (16)  is 

0^S,<:<M.  (18) 

Also,  define  a (m  + 1)  x (m  + 1)  block  matrix  , 


O 

0 

1 

o 


o 


'T  o 

0 I 


(19) 


where  we  require  m ^ p ^ 1 for  this  definition.  Then,  using  the  nota- 
tion established  in  equation  (18), 


(20) 


where 
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If  -- 


At,  l=  + ^F 


Also, 


0 ^ sej 

|Qt’"]  = iir,  -for  0£r,s<rn- 

Then,  the  r,u-th  block  of  the  product 

lQr"corL--^j<?rLf«r),iGrL 

ts„i-s„ArK':&..r-s..Ar] 

t\}„  \X..  - S-.  Ar" 

+j.s.A:x'':Ar"] 


(21) 


(22) 


5;f*  0 


(23) 

=xvs„^Ar«t:-  Ar%xj»  r 
■ k::  - s,.  - s.  i k*'!  k s«  cc 


In  the  last  line,  above,  we  have  used  equation  (21)  to  simplify  equation  (.23). 
At  this  point,  we  consider  four  subcases: 

(a)  for  1 ^ r,  u ^ m,  equation  (23)  reduces  to 

(b)  for  r = 0,  u = 0,  equation  (23)  becomes 
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(24) 


but  by  use  of  equation  (14)  , the  sum  on  s in  the  last  term  of  equation 
(24)  is  , in  which  case  the  last  two  terms  of  equation  (24)  cancel. 
We  are  left  with 


using  equations  (12)  and  (95)'; 

(c)  for  r = 0,  1 < u ^ m,  equation  (23)  yields 

±Arc:  -0. 

using  equation  (14);  and 

(d)  for  u = 0,  1 ^ r ^ m,  equation  (23)  yields 

-r  ^ » I 


(25) 


(26) 


(27) 


since  this  is  the  conjugate  transpose  of  equation  (26).  Therefore,  we 
have 


Up  0 o 

0 1^: 

0 K'f 


0 


0 'R 


(p) 


K 


0 


0-  0 


(28) 
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This  relation  holds  for  m 2;^  p ^ 1 , as  noted  under  equation  (19)  (some 
relations  for  determinants  are  noted  in  appendix  B) . 

Now,  let  be  arbitrary  nonzero  complex  M x 1 column  matrices. 

Then,  using  equation  (28), 


I 

v; 

1 

1 

y-. 

(29) 


We  recall  that  Up  is  positive  definite,  by  the  previous  section. 
Therefore,  if  |R^P|  is  positive  definite,  then  ^m^^  posi- 
tive definite,  which,  in  turn,  implies  that  is  positive  definite. 

That  is,  for  m ^ p ^ 1, 


if  ig  positive  definite,  then  <R^P^  is  positive  definite . (30) 

In  particular,  letting  m = p,  we  see  that  if  is  positive  defi- 
nite, then  (R^P^  is  positive  definite.  But  by  equation 

(17).  Hence,  if  is  positive  definite,  then  ^^P)  is  positive 

definite.  But  «(0?  = Rq  is  positive  definite  (see  equation  (15)). 
Therefore,  we  conclude  by  induction  that 


«(P) 

P 


is  positive  definite  for  all  p. 


(31) 


This  statement  is  used  as  a priori  information  in  Burg's  derivation  in 
the  known  correlation  case  (see  reference  3,  page  85) . 


Now,  we  return  to  equation  (30)  with  this  information  and  can  draw 
the  conclusion  that  ^CP)  is  positive  definite  for  all  m ^ p.  Finally, 
using  equation  (17) , 0e  can  state 

is  positive  definite  for  all  m and  p.  (32) 

For  fixed  p,  since  is  positive  definite  for  all  m,  (the  ele- 

ments of)  f^P^  cannot  tend  to  infinity  as  m tends  to  infinity,  since 
R^P'  = Rq  is  fixed.  Therefore,  recursion  (11)  is  stable.  This  implies 
(using  equation  (23) 1)  that 
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a,+(x-  - w %%) 


(33) 


possesses  all  its  zeros  inside  the  unit  circle  in  the  z-plane;  that  is, 
predictive  error  filter  (z)  is  minimum  phase. 

The  proof  above  hinges  critically  on  the  positive  definiteness  of 
U , which  was  demonstrated  in  the  previous  section.  In  particular,  this 
condition  is  employed  in  equation  (29)  to  guarantee  that  the  right-hand 
side  be  positive. 

A word  of  caution  about  an  apparent  alternative  proof  is  worth 
mentioning  here.  Having  shown  that  Up  is  positive  definite,  one  might 
be  tempted  to  define  r(P^  by  the  inveree  of  equation  (165)  ,1 


(34) 


according  to 


J. 

in 


= I ex|.(i  2-rfw,a)  (f)  , all 


(35) 


24 


It  is  obvious  that  G^P^(f)  in  equation  (34)  is  positive  definite  for 

any  f;  and  it  is  now  easy  to  demonstrate  that  jJ(P)  ij 

m 


Is  positive  definite: 


[r 


% 


^ V. 


5;T*0 


J. 


Ik 


(36) 


-L 


" 24 


since  G^P^(f)  is  positive  definite  for  any  f, 
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I 

i 

) 

{ 

However,  the  problem  is  that  we  now  would  have  to  show  that  | 

as  generated  by  equation  (35),  satisfies  the  recurrence  (11).  An  example  j 

in  appendix  C shows  that  for  an  unstable  sequence,  the  values  returned  j 

by  equation  (35)  are  not  the  same  sequence;  thus,  equation  (35)  should  | 

not  be  used  until  after  the  stability  of  (R^PM  has  been  ascertained. 

1 

ALIASED  CORRELATIONS  VIA  FFT  1 


Based  upon  the  previous  results,  we  know  that  we  can  express 

^ A^.XeAc|>(-,'2Trf*ria)K^  , Ifl  ^ il  ^ 


and 

_L 

- J df  Gr ("f ),  all 

~ zk 

We  have  dropped  the  superscript  p above,  since  the  results  to  follow 
will  hold  for  any  correlation-spectrum  pair  satisfying  equations  (37) 
and  (38) . 

If  spectrtun  G(f)  is  calculated  only  at  a discrete  set  of  N + 1 
/ 1 1 \ 

points  on  I-  (vdiich  is  a typical  practical  situation  for  plot- 

ting purposes,  for  example),  a discrete  approximation  is  afforded  to 
the  integral  in  equation  (38).  It  is,  for  trapezoidal  weights  (wj, } , 

That  is,  the  discrete  approximation  to  integral  (38)  yields  aliased  sam- 
ples of  correlation  sequence  {Rn,}  at  separations  of  Np;  this  is  easily 
proven  by  substituting  equation  (37)  into  the  left-hand  side  of  equation 
(39)  and  interchanging  summations. 

» 

The  aliased  sequence  {R^i)  has  period  Np.  Therefore,  is  a good 

approximation  to  R^  for  |m|  < Np/2  if  |R^|  is  sufficiently  small  for 


(37) 


(38) 
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|m|  > Np/2.  (Generally,  Np  >>  P3£5i-  in  the  linear  predictive  approach, 
and  this  is  true.)  The  reason  for  considering  this  approach  to  the 
approximate  evaluation  of  correlation  sequence  (R,„}  follows. 


The  left-hand  side  of  equation  (39)  can  be  accomplished  by  means  of 
an  Np-point  FFT  (one  FFT  for  each  element  of  the  M x m matrices  involved). 

For  trapezoidal  weights,  using  the  fact  that  cf-  equation 

(39)  is  expressible  as  ' ' 


»r-Np/2 


[i  2Trhvn/^f'') 


©K|>(|2irU/Np)G; 


(40) 


where  we  have  defined 


Gr^  = ^(iv)  > ^ 


(41) 


Letting  n = Np  + m in  the  first  sum  of  equation  (40),  and  n = m in 
the  second  sum,  we  obtain 


tXj>(i2ir«iw/Np)X,  , 


(42) 


where  M * m matrix 


1 •- 


t z 


(43) 


But  equation  (42)  is  recognized  as  an  Np-point  FFT  of  the  matrices 
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^o,  - > ^ , •••/  ^-1  ' 


(44) 


thus,  we  obtain  Rq,  Ri,  . • R(g^_|  by  means  of  tiis  Np-point  FFT,  one 

FFT  for  each  element  of  the  M x m matrices.  (The  quantities  (Rni^ 

|m|  < Np/2  are  available  by  use  of  the  periodic  nature  of  sequence 
(Rm^^  of  ^ Np-point  FFT  to  obtain  (good)  estimates  of 

correlation  sequence  {R  } circumvents  the  use  of  recursion  (11),  which 
would  yield  the  exact  correlation  sequence  {R,j}  . It  can  save  time  in 
some  cases  and  uses  already  available  quantities  {Gj^}  , if  they  have 
been  computed  previously  for  plotting  or  observation  purposes. 

REAL  PROCESSES 

The  preceding  results  for  complex  multivariate  processes  can  be 
specialized  to  real  processes.  We  have,  from  equations  (171)^  and  (39), 


\ 


'K  > 'H 

Therefore,  equation  (39)  becomes 


1. 


(45) 


K,  H '■»r  ('  fij, , 


Nf' 


(46) 


where 


, y-O  ov  Mf/a 
0-  k < N/p/z 

Now,  let  the  elements  of  matrices  Gj^  and  R^  be  expressed  as 

(tt) 

Then,  is  real  for  all  t;  and  from  equation  (46), 


(47) 


(48) 


'R' 


(w) 


Nf/a 


Co5[2-irk»i«//\^:)Grj^ 


Np  A K = o 


(49) 
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In  addition,  since 

a!:’,  ^ 


^ (itf) 


itt) 


(50) 


n(^^) 


the  fundamental  range  of  m is  [0,  Np/2]  for  sequence  }. 

REAL  BIVARIATE  PROCESSES 

We  can  specialize  further  to  the  bivariate  case,  M = 2,  and  make 
use  of  some  of  the  properties  previously  discussed.  (The  goal  of  these 
manipulations  will  not  be  clear  until  the  final  result.)  Define  the 
complex  scalar  sequence  {Uj^}  such  that 


K 


(51) 


Then , 


hV-l 

X 

k-  0 

\ 


\ ex  ^ (^1  i 2Trk*»i/Np^ 

[GrJ%  i ( 2irk*w/N,) 


k--o 

^ +i  \ 2iry,wlH^ 


(52) 


If,  on  the  right-hand  side  of  equation  (52),  we  let  n = k in  the 
first  sum,  snd  n ~ Np  - k in  the  second  sum,  v/e  get 


13 


the  last  step  by  equation  C49) ; that  is,  using  equation  (52)  again, 

A-'  r 


(54) 


Thus,  one  Np-point  FFT  of  scalar  sequence  defined  in  equation 

(51).  will  give  both  (aliased)  real  scalar  autocorrelations  and 

and  by  the  statement  under  equation  (50), { } need  be 

printed  out  only  for  0 ^ m ^ Np/2. 


For  the  crosscorrelation,  equation  (46)  yields 
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Ks  0 

H./2 


(55) 


,'\N,/2 


This  Np-point  FFT  of  ~ + 1 nonzero  numbers  would  yield 
and  from  equation  (39) , since 


Np-l 


R = R (for  general  complex  M x M matrices), 
-m  m 


(56) 


it  follows  (using  the  periodicity  of{R^  that  for  the  present  case 

(57) 
^F 

Thus,  print  out  of  and  for  0 £ m ^ ^ suffices  to  give 

complete  information  about  the  aliased  crosscorrelation.  Furthermore, 
all  this  information  is  available  from  the  single  Np-point  FFT  of  equa- 
tion (55) . 


In  summary,  only  the  two  FFTk  indicated  in  equations  (54)  and  (55) 
need  be  conducted  to  obtain  complete  information  about  the  aliased 
correlation  sequence  (fin,},  for  M = 2.  These  relations,  in  addition  to 
the  exact  correlation  recursion  (11),  have  been  incorporated  in  the  1 
FORTRAN  program  listed  in  appendix  D.  The  comments  in  appendix  K of 
the  earlier  report^  are  relevant  here  also. 
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SU^*^ARY 


It  has  been  shown  above  that,  for  the  weighting  introduced  in 
equation  (136),1 

= Vp-,  ^ cV)oice  2,  (58) 

Up  and  Vp  are  guaranteed  positive  definite,  and  the  correlation  recur- 
sion (11)  is  stable.  Therefore,  equation  (58)  is  a sufficient  condition 
for  the  desired  properties  to  hold  true.  It  is  not  known  whether  this  is 
a necessary  condition,  that  is,  whether  equation  (58)  is  the  only  choice 
that  results  in  the  desired  properties  of  positive  definiteness  and 
stability. 


However,  for  M = 1,  since,  by  equation  (129),^  1)  = Vp_j , it  is 

possible  to  show  that  ^ 


Ap.,  * ^p-,  0 


(59) 


is  the  only  choice  that  guarantees  the  desired  properties  (see  refer- 
ence 1,  page  32).  Namely,  equations  (124), ^ (130), 1 and  (114)1  yield 
scalar 


In  addition,  if  the  data  samples  happen  to  take  on  values  such  that 


(60) 


(61) 


then 


(62) 


the  sample  mean  of  the  original  data  is  (made)  zero,  this  choice 
is  not  possible  for  p = 1.  For  p > 1,  the  sample  means  of  and 

are  not  necessarily  zero. 
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which  is  always  larger  than  1 (unless  Ap.j  = Tp.j) ; then^Up  is  negative 
and  an  unstable  correlation  recursion  results.  tTius,  equation  (59)  is 
the  only  choice  that  guarantees  positive  Up  and  a stable  correlation 
recursion,  regardless  of  the  data  set,  for  M = 1. 

It  should  be  noticed  that  the  absolute  level  of  the  weights  is  not 
specified  by  equation  (59).  Thus,  for  M ^ 2,  freedom  in  equation  (58), 
at  least  to  the  extent  of  a common  scale  factor,  must  he  allowed. 
Whether  this  is  the  only  degree  of  freedom  allowed  to  the  choice  of 
A , and  r , is  unknown  for  M > 2. 

p-1  p-1  - 


17/18 
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Appendix  A 


SOME  PROPERTIES  OF  COMPLEX  MATRICES 


An  arbitrary  complex  square  matrix  A is  called  real  definite 

= r C'^ol)  W Qwy  % 

where  V is  a complex  column  matrix. 


if 

(A-1) 


It  then  follows  that 


A real  definite  A^^  = A,  real. 


(A-2) 


where  are  the  eigenvalues  of  A. 


For  proof,  first  take  the  conjugate  transpose  of  equation  (A- 

V**  h ^ ^ ^ ■ 

Subtracting  equations  (A-1)  and  (A-3)  gives 

~ 0 W anj  ^ 

Therefore, 

A^-A=  0,  or  A'’  - A 

Also,  if  {V]^}  are  the  eigenvectors  of  A,  then 

AV,  = \Vk, 

v;ak  = \v:v,  . 


1). 

(A-3) 


(A-4) 


(A-5) 


(A-6) 


Since  the  left-hand  side  and  are  real,  Xj^  is  real, 


If  r in  equation  (A-1)  is  positive  for  any 'lMo,then  A is  said 
positive  definite.  It  follows  that 


to  be 


A positive  definite  A 


H 


A,  {X^}  > 0. 


(A-7) 


A-1 


f 


\ 

I 


TR  5729 


The  proof  is  the  same  as  the  proof  above,  except  that  now  V^AV,  > 0 in 
equation  (A-6) . 

Now,  we  are  in  position  to  prove  that 

A positive  definite  1 

B positive  definite  J ate  nil  positive. 

For  proof,  let  {Xj^}  and  {Vj^}  be  the  eigenvalues  and  eigenvectors  of  AB; 
then,  we  have 

- \kV, 

^x=\A''v. 

v:bx.\v:a'\=\.i;\'’vo'’a(a'v.\ 


where  we  have  used  A"  = A (equation  (A-7)).  Since  A and  B are  positive 
definite,  the  left-hand  side  and  the  factor  multiplying  Xj^  are  positive. 
Therefore,  Xj^  is  positive. 

It  should  be  noted  that  AB  need  not  be  Hermitian  or  positive  defi- 
nite. For  example,  if 


o<  0*']  j , 

- I •<  oc>  0,  ) 

-1 

M"'  u real,  ^>0,  h > 

)*  J ^ 


(A-10) 


then. 


o<y^3 

AB  U (A-11) 

Since  the  main  diagonal  terms  of  AB  need  not  be  real,  AB  is  not 
necessarily  Hermitian.  Also,  if  we  assume  that  AB  is  positive  definite, 
equation  (A-7)  says  that  AB  is  Hermitian,  which  is  contradictory. 


A numerical  example  follows: 
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-2  i-r 

B = 

'i 

l+-i' 

li-i  2 , 

i-i 

2 . 

(A-12) 


A and  B are  positive  definite  and  Hermitian.  The  eigenvalues  of  both 
are  {Xj^}  = 0.  Their  product  is 


4- >2  4 

. 4 4+i2j  ’ 

with  eigenvalues  4 ± 2 .V?  > 0,  as  predicted, 
nor  positive  definite  because,  for  instance, 


(A-13) 


But  AB  is  not  Hermitian 


fi 


4-i  2. 


(A-14) 


The  matrix  AB  in  equation  (A-13)  points  out  that  specifying  a 
matrix  to  have  positive  eigenvalues  does  not  make  that  matrix  positive 
definite.  However,  if  the  matrix  is  also  Hermitian,  we  have  the  genera- 
lization of  equation  (A-7)  to 

A positive  definite  ^ A^  = A,  {Xj^}  > 0.  (A-15) 


A-3/A-4 
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Appendix  B 


RFLATinNS  OF  DETERMINANTS 


Since  det  = 1 (see  equation  (19)),  equation  (28)  yields 

^ dct  0^  (RJ.,  , '*1  - p 


(R-1) 


Setting  m = p in  equation  (B-1)  and  employing  equation  (17),  there  fol- 
lows 


(R^’  - det  U,  del  . 


Since  = Rq  = IJq  (see  equation  (95)^),  this  recursion  may  be 

written  in  closed  form  as 


det  (Rp  = det  V)^ 


(B-3) 


This  relation  is  given  in  Burg,^  page  86. 


By  letting  m = p ♦ 1,  p + 2,.  . ..in  equation  (B-1),  it  follows 
immediately  that 

dtt  (RJI  = (det  V)^)  ^ det  , Vn  > p. 


In  addition,  for  m < p,  using  equations  (17)  and  (B-3), 

dettRl"  . detC-  itdetU,  »<  r 

K=o 

Combining  equations  (B-4)  and  (B-5) , we  have 


(B-5) 


det  e / 


I t^  de+  , m<  p 

' K=o 

“iVdet  U 


B-l/B-2 
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EXAMPLE  OF  UNSTABLE  CORRELATION  RECURSION 


Consider  the  univariate  (M  = 1)  correlation  values, 
R = r^"’L  all  m,  r real  and  positive. 


(c-n 


The  value  of  r can  be  greater  or  less  than  unity.  The  z-transform  of 
equation  (C-1)  is 


^ 'h  hi»-t 


(C-2) 


Now . 


5,  = 


-2 


s = 


i 2' 


>T 

iT  1^1 


(C-3) 


But,  if  r > 1,  there  is  no  common  region  of  convergence;  also, 
sequence’cR^}  i?  unstable  if  r > 1.  Nevertheless,  if  we  blithely  add 
terms  in  equation  (C-2),  we  get 


(C-4) 


Then,  continuing  on,  setting  z = exp(i2TTfA)  and  multiplying  by  A 
[©Xf(i  ?Tt  f a) -1^  (i  “ r] 


(C-S) 


which  is  real,  and 


C-1 


TR  5729 


± 


^2*"  ^ ^ ~ ^ 


_x 

u 


Iklii^ 

C«>«|« 


,.,  * (l-rX^-i-) 


In  the  following,  let  r ^ 1,  a = min(r,i),  and  6 = max(r,^). 


Jwl 


■fov'  qH  hi 


This  is  a stable  sequence  for  any  r.  But,  notice  that  if 

r<  1 ; ‘^  = ^ ~ "r  ^ ~ **•» 


whereas,  if 


(C-6) 

Then, 

(C-7) 

(C-8) 

(C-9) 


The  former  sequence  is  correct;  the  latter  is  not.  Yet  both  are 
stable.  So,  although  equation  (C-6)  always  generates  a stable  sequence, 
it  is  not  necessarily  the  original  sequence. 
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Appendix  D 

FORTRAN  PROGRAM  FOR  SPECTRAL  ANALYSIS 


A FORTRAN  listing  of  the  spectral  analysis  technique  is  given  in 
this  appendix,  in  addition  to  a sample  printout  of  ^ application.  The 
notation  and  scaling  adopted  is  identical  to  that  given  in  reference  1, 
appendix  K.  The  equation  numbers  referenced  are  those  in  the  earlier 
report,!  except  in  Subroutine  ACM,  where  they  correspond  to  the  equa- 
tions in  this  report. 
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$kiiA(M»M)  »WC(f  »M)  ,WD(^1,y)  »WE(V»M) » AIC  (PMAX)  , AlCO  (2)  »S(M) 

EQUIVALENCE  ( X » Y ) » ( /^  IC  ( 1 ) * AICO  ( 2 ) ) 


FORMAT(I10.E20,8) 
PRINT  12 t PBEST 
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